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Abstract. We present the final results obtained by the MACRO experiment in the search for GUT magnetic 

monopolos in the penetrating cosmic radiation, for the range 4 x 10~' <(]<!. Several searches with all the 
MACRO sub-detectors (i.e. scintillation counters, limited streamer tubes and nuclear track detectors) were 
performed, both in stand alone and combined ways. No candidates were detected and a 90% Confidence 
Level (C.L.) upper limit to the local magnetic monopole fiux was set at the level of 1.4 x 10~^® cm~^ s~^ 
sr~^. This result is the first experimental limit obtained in direct searches which is well below the Parker 
bound in the whole (3 range in which GUT magnetic monopoles are expected. 
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1 Introduction 

Within the framework of Grand Unified Theories (GUT) 
of the strong and electroweak interactions, supermassive 
magnetic monopoles (MMs) with masses m ^ 10^^ GeV/c^ 
could have been produced in the early Universe as intrin- 
sically stable topological defects at the phase transition 
in which a simple gauge symmetry left an unbroken U(l) 
group 0]. At our time they can be searched for in the 
penetrating cosmic radiation as "fossil" remnants of that 
transition. The detection of such a particle would be one 
of the most spectacular confirmation of GUT predictions. 

The velocity range in which GUT magnetic monopo- 
les should be sought spreads over several decades. If suffi- 
ciently heavy (m > lO^^GeV/c^), GUT magnetic mono- 
poles would be gravitationally bound to the galaxy with 
a velocity distribution peaked at /3 = v/c ~ 10""^ 
MMs trapped around the Earth or the Sun are expected 
to travel with P ~ 10^^ and ~ 10~^, respectively. 

Intermediate mass MMs (10^5 < m < lO^"^) could 
have been produced in later phase transitions in the early 
Universe at a lower energy scale lighter magnetic 
monopoles, with masses around 10''-;- 10^*^ GeV/c^, would 
be accelerated to relativistic velocities in one or more 
coherent domains of the galactic magnetic field, or in the 
intergalactic field, or in several astrophysical sites like a 
neutron star 

Theory does not provide definite predictions on the 
magnetic monopole abundance. However, by requiring 
that MMs do not short-circuit the galactic magnetic field 
faster than the dynamo mechanism can regenerate it, a 
flux upper limit can be obtained. This is the so-called 
Parker bound (~ 10~^^ cm~^s~^sr~^ Q), whose value 
sets the scale of the detector exposure for MMs search. 
The original Parker limit was re-examined to take into 
account the almost chaotic nature of the galactic magnetic 
field, with domain lengths of about 1 10 kpc; the limit 
become mass-dependent. An Extended Parker Bound 
(EPB) at the level of 1.2 x W-'^^{m/10^'^) cm-^s-^sr-i 
was obtained p^ . 

MACRO was a large multipurpose underground de- 
tector located in the Hall B of the Laboratori Nazionali 
del Gran Sasso (Italy); it was optimized for the search 
for GUT magnetic monopoles with velocity /? > 4 x 10"^ 
and with a sensitivity well below the Parker bound. The 
detector, which took data up to December 2000, is fully 




Fig. 1. A cross-sectional view of the MACRO detector showing 
the layout of the 3 sub-detectors and of the rock absorber. 



described in |5[g]. It was arranged in a modular structure, 
it was divided into six 12.6 x 12 x 9.3 m"^ sections referred to 
as SuperModules (SM), each one with separate mechanical 
structure and electronics readout. The detector's global di- 
mensions were 76.5 x 12 x 9.3 m^ and its total acceptance 
for an isotropic flux of particles was ~ 10,000m^sr. 

Redundancy and complementarity were the primary 
features in designing the experiment. Since we could not 
reasonably expect more than a few MMs during the 
detector lifetime, we deemed crucial to have multiple 
signatures and the ability to perform cross checks among 
various parts of the apparatus. To accomplish this, the 
detector consisted of three independent sub-detectors: 
liquid scintillation counters, limited streamer tubes and 
nuclear track detectors, each of them with dedicated 
and independent hardware. In Fig. ^ a cross sectional 
view of the apparatus is shown. Also visible are the 
seven horizontal absorber layers (which set at ~ 1 GeV 
the minimum energy threshold for throughgoing muons); 
notice the separation of the detector into a lower and an 
upper detector (the attico). 

The background that magnetic monopole searches 
have to fight with is mainly due to muons of the cosmic 
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radiation and natural radioactivity. Thus large detectors 
have to be installed in underground laboratories. At Gran 
Sasso the minimum thickness of the rock overburden 
above the detector is 3150hg/cm^. The cosmic radiation 
muon flux in Hall B is ~ lm~^h~^, almost a factor 
10^ smaller than that at the surface (only muons with 
minimum energy of ~ 1.3 TeV can cross the mountain and 
reach MACRO) with an average muon residual energy of 
~ 300 GeV §. 

The signatures of the passage of a GUT magnetic 
monopole across the detector depend strongly on its 
velocity For this reason, different hardware systems 
were designed and operated to give optimum sensitivity in 
different /3 values. Different analysis strategies were also 
adopted, depending on the (3 range of interest and the 
subdetector(s) used; the entire MM fi range (4 x 10~^ < 
/? < 1) was thus covered. In the different analyses we took 
into account both the MM signatures and the background 
characteristics. 

The uniqueproperty of a fast magnetic monopole with 
(3 > 10~^ is its large ionization power compared 

either to the considerably slower magnetic monopoles or to 
the minimum ionizing electrically charged particles. The 
searches for fast MMs look for large energy releases in the 
detectors; the background is mainly due to high energy 
muons with or without an accompanying electromagnetic 
shower. On the other hand, slow magnetic monopoles 
should leave small signals spread over a large time window; 
a /3 ~ 10~^ monopole could have a Time of Flight (ToF) 
across the detector as large as ~ 1 ms. This implies the use 
of specific analysis procedures that allow the rejection of 
the background mainly due to radioactivity induced hits 
and possible electronic noise. 

In this paper we report the final results of several ma- 
gnetic monopole searches performed with MACRO (Sect. ||) 
some early results were already published in Q. In each 
Section we discuss the analysis criteria together with the 
results of each search; technical details may be found in 
various papers fully describing the procedures and their 
application to the first data samples [|l0yil|Jl^ , |l3| , p^ . In 
Sect. H the result of the combination of all the various 
searches is reported. In order to compare the MACRO re- 
sults to those of other experiments or to theoretical mod- 
els, we present upper limits for an isotropic flux of bare 
MMs with magnetic charge equal to one Dirac charge 
9 = go = e/2a, (where e is the electron charge and a 
the fine structure constant) and nucleon decay catalysis 
cross section smaller than 1 mb . These aspects are dis- 
cussed in Sect. |[ the conclusions are given in Sect. I A 
dedicate analysis for the search for magnetic monopoles 
accompanied by one or more nucleon decays along their 
path was also performed. The results are reported in a 
separate paper fSTt]. 



2 Experimental searches 
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Fig. 2. Results of the measurements of the Slow Monopole 
Trigger sensitivity (black squares) and of the WFD analysis 
efficiency (stars) plotted vs. /3. The connecting lines serve only 
as a guide to the eye. The expected light yield of a MM is also 
plotted. 



2.1 Searches with Scintillators 

The MACRO liquid scintillator system was organized in 
three layers of horizontal and four layers of vertical coun- 
ters, as shown in Fig. |l|. For atmospheric muons cross- 
ing the apparatus this system provided particle position, 
energy deposition and ToF resolutions of about 11 cm, 
1 MeV and 700 ps, respectively. The response of liquid 
scintillators to heavily ionizing particles was studied both 
experimentally |15| and theoretically p6| , p7| , p^ and their 



sensitivity to particles with (3 down to 
measured. 



10 ^ras directly 



2.1.1 Wave Form Digitizer (WFD) analysis 

Slow moving magnetic monopoles (in the range 10^^ (3 "Z 
4.1 X 10""^) were searched using dedicated hardware, the 
Slow Monopole Trigger (SMT) and a 200 MHz custom- 
made Wave Form Digitizer (WFD) system. The main goal 
of the SMT was to remain sensitive over the entire range 
of widths and amplitudes of pulses (which could as well 
be just a train of single photoelectron pulses lasting over 
several microseconds) expected for slow moving magne- 
tic monopoles while suppressing efficiently narrow (10 - 
50 ns) pulses due to isolated radioactivity and cosmic 
ray muons. The SMT by itself offered only a hit regis- 
ter that recorded the scintillator counters that satisfied 
the trigger conditions. The photomultiplier tube (PMT) 
pulse shape information was recorded by a custom made 
200 MHz WFD system. Redundant time-of-flight (TOF) 
information for every candidate event was recorded by a 
stand-alone TDC. The design and implementation of the 
triggering and recording electronics used by this search 
has been described in detail elsewhere . 

The heart of this slow MM search was the analysis 
of the digitized PMT signals as recorded by the 200MHz 
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WFDs. Several wave form analysis methods have been 
employed in the search for slow MM searches in the initial 
data taking and results from these have been reported 
elsewhere [0 . In this final search the SMT hardware |^,^ 
was used as the pattern recognition method for identifying 
a magnetic monopole candidate event [ p^ . 

The sensitivity of the trigger hardware as well as of the 
entire slow MM analysis was measured experimentally in 
situ and multiple times throughout MACRO'S live time 
via LED light injection in each of the scintillation coun- 
ters. The LED calibration system |^,D allowed the gener- 
ation of a magnetic monopole-like signal in the detector of 
arbitrary direction, ionization yield and velocity. A grid of 
LED pulse widths and heights allowed the generation of 
MM-like PMT pulses on a channel-by-channel basis. The 
corresponding /3-range is 10"'' < (3 < 5 x 10~^ and the 
light yield from few photoelectrons up to several times the 
yield of a muon. The generated wave forms were recorded 
by the WFD system independently of whether the SMT 
fired or not. Then, an off-line wave form analysis (the same 
used to analyze the real events) established the ratio of 
the events that the SMT had selected over the ones ex- 
pected. This yielded the trigger sensitivity as a function 
of the particle's velocity and it is shown in Fig. |^. More 
than 95% of the detector's channels exhibited efficiency 
above 99% for light yields greater or equal the ones de- 
fined by the curve. The generated MM-like signals were 
then fed into the analysis wave form selection algorithm 
in discussion. The stars in Fig. ^ show the analysis ef- 
ficiency, when more that 99% of the simulated MM-like 
wave forms -corresponding to lights yields and velocities 
across the expected range- were successfully identified by 
the analysis algorithm. 

There were three main instrumental sources which cau- 
sed the hardware trigger to be less efficient than a software 
one: (a) electronics were operating at a very low threshold 
(2 mV) making them susceptible to ground loops and 
low amplitude interference; (b) electronics suffered from 
intrinsic pulse streching that may overestimated pulse 
duration; (c) helium and hydrogen contamination in the 
PMT envelope increased its activity, leading sometime 
to overestimate of pulse duration. A software pattern 
recognition method corrects for them, and the efficiency 
of this algorithm is folded in the trigger efficiency already 
discussed. 

Magnetic monopole candidate events were selected 
as follows. The hardware trigger (SMT) was required 
to be present in at least two of the detector faces and 
within 1msec. This rejects mainly noise and radioactivity 
induced events. The two layers requirement reduces by 
17% the detector's acceptance. The TOF requirement of 
1 msec was dictated by the hardware, namely the mode 
of operation of the WFD system [gjH]. The presence of 
the scintillator's fast particle trigger was then used as a 
veto. The fast coincidence (1 /is) between faces required 
by this trigger selects particle tracks corresponding to /3 
well above the SMT's sensitivity (Fig. ^). It rejects mainly 
cosmic ray muons that triggered the SMT either due to 
pulse stretching or afterpulsing. A final loose cut on the 



hits per face (less than 4) and number of layers (less 
than 5) was applied to reject events induced by electronic 
noise, multi-muon or muon shower events. None of these 
requirements is expected to affect the sensitivity to slow 
MMs while any effect on the acceptance has been taken 
into account in the Monte Carlo. After applying the above 
cuts on the dataset of ~4.75 years of data-taking and ~28 
million SMT's, there were 35901 events left. These were 
attributed to cosmic ray muons (« 1/3) and electronic 
noise (~ 2/3) leaking out of our fast particle and loose 
noise/high multiplicity veto. These events were further 
analyzed by our wave form analysis method. 

The wave form analysis required that PMT signals 
satisfied the pattern recognition criteria on an end-to- 
end, face-to-face basis. 40 events survived. (There was an 
initial running of the detector that suffered from wave 
form memory saturation and effectively wave form loss. 14 
events involving primarily large ionizations were allowed 
to pass this stage of the analysis). 

The 40 remaining candidates were individually ana- 
lyzed using all available information from the scintillator 
sub-system (the wave form shape, timing information from 
the wave forms, timing information from the precise muon 
system -ERP- and from the SMT TOF system) None 
of the events was compatible with the passage of a slow 
MM. Instead, based on information available from both 
the scintillator and the streamer tube sub-systems, the 40 
events were classified as follows: 

— 25 as bipolar electronic noise; using primarily the WFD 
information, they are distinguished from a genuine MM 
signal due to the preciable amount of positive content 
in the waveforms (PMTs operated with negative volt- 
age). 

— 10 as muon-induced showers; using primarily the fast 
muon and occasionally WFD TOF information, these 
events reflect fast crossing of the apparatus that failed 
to fire the fast particle veto. 

— 4 as probably spurious electrostatic discharges of the 
LED system; using primarily geometrical arguments 
refiecting the inconsistency of the event with the 
passage of a monopole crossing the detector at a 
straight line at a constant velocity, and finally, 

— the "spokesmen even" . This event belongs to a par- 
ticular run in which a LED-induced event (simulat- 
ing a close to a realistic f3 ~ lO^'^ MM) was inten- 
tionally and secretly generated to test the efficiency of 
our analyses to detect a MM signal. The two analyses 
that should have been sensitive to the faked "spokes- 
men" MM (the present analysis, and the PHRASE 
one), found the event. The "spokesmen event" was a 
LED-generated event and was subject to the limita- 
tions the LED system had in generating a MM-like 
pulse, namely its slow rise-time. 

No event was compatible with the passage of a single slow 
magnetic monopole. The analyzed data were collected 
in the period July 1995 - May 2000, corresponding to 
~ 4.3 live-years. The detector performance during the 
data-taking period was monitored on a run-by-run basis. 
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Fig. 3. Calculated crossing time in a scintillation counter 
versus pulse width for PHRASE monopole candidates in the 
range 1.2 x 10"^ < /3 < 5 x 10"^ (for clarity reasons only 40% 
of the candidates is shown). The diagonal line indicates the 
minimum pulse width expected (corresponding to a MM with 
/3 ~ 5 X 10~^) for 15 cm path length. No candidates meet or 
exceed this expectation. 



The total exposure was / Adt = 0.9 x 10^^ cm^ssr. By 
taking into account the trigger and analysis efBciencies 
we obtained a 90% C.L. flux upper limit of 2.5 x 
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"WFD" in Fig. 



for 10-^ 

i9|. 



< /3 < 4.1 X 10' 
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2.1.2 The PHRASE analysis 

Magnetic monopoles in the velocity range 1.2 x 10" 
f3 < 10~^ were searched for using the PHRASE (Pulse 
Height Recorder And Synchronous Encoder) system. It 
was a trigger and energy reconstruction processor which 
generated a trigger condition when the energy deposition 
in a scintillator was larger than a preset threshold, Eth ~ 
7MeV [|j20|2l|. 

The energy deposition was reconstructed by integrat- 
ing the photomultiplier pulses (recorded by a 100 MHz 
Wave Form digitizer system), taking into account their 
relative amplitudes, the response functions of the scintil- 
lation counters and the liquid scintillator, photomultiplier 
and electronic saturation effects. All these points were 
studied in detail by using laser light of various intensi- 
ties and atmospheric muon pulses; appropriate algorithms 
were developed for correcting any nonlinearity which may 
be present in large energy loss rates as those expected 
from magnetic monopoles in the analyzed grange. The 
analysis technique is described in detail in pJ^ . 

The magnetic monopole candidate selection required 
hits in two scintillator layers and a minimum energy 
deposit of 10 MeV in each of them, spread out over 
no more than four adjacent boxes (this cut removed 
most of the natural radioactivity background but not 
atmospheric muons). The 100 MHz sampling frequency 



of the PHRASE Wave Form digitizers and the 10 MeV 
energy cut defined the minimum MM velocity which this 
search was sensitive to: /3min = 1-2 x 10~^. The velocity 
of a particle was reconstructed by using the TOF between 
the crossed scintillator layers; a minimum distance of 2 m 
between the hits was required to ensure accurate timing 
reconstructions. The timing and position uncertainties 
produced a tail in the muon velocity distribution which 
is peaked aA, (3 — 1. Atmospheric muons were rejected 
by setting the (3 upper limit for the magnetic monopole 
analysis to 0.1. Finally, a minimum path length of 15 cm 
was required in each scintillation counter (the typical path 
length for a crossing cosmic ray muon was > 20 cm). All 
geometrical cuts were taking into account in computing 
the detector acceptance. 

The candidates surviving the selection (a few thou- 
sands in ~ 11 years of running) were grouped in two 
(3 overlapping sub-ranges, medium (1.2 x 10~^ < (3 < 
5 X 10^3) and high (4 x lO'^ < (3 < lO'^) velocities. This 
was motivated by the fact that the MM light yield and 
crossing time (with respect to muons) differ significantly 
in these two regimes. Crossing muons released ^ 40 MeV 
in each traversed scintillation layer, with a pulse width of 
~ 35 ns due to the convolution between the time profile 
of the scintillation light emission and propagation and the 
PMT time jitter. The pulse width expected from a magne- 
tic monopole in the medium P range, corresponding to the 
time needed to cross a counter, is 100-^400 ns. In the high 
(3 range the expected counter crossing time is 6 120 ns; 
thus the corresponding pulse width is comparable to that 
of a muon. The expected energy loss rate, for MMs in the 
medium /3 range, is 10 30 times larger than that of a 
muon, and 30 ~ 60 times for MMs in the higher f3 range. 

For candidates with 1.2 x 10~^ < /3 < 5 x 10~^ we 
compared the measured pulse width with the expected 
counter crossing time. The counter crossing time was 
calculated from the measured /?, assuming a minimum 
pathlength in the scintillation counter of 15 cm (Fig. ||). 

For a particle crossing the detector with the measured 
velocity these two numbers must be consistent. No final 
candidates satisfied this condition with the exception 
of the LED simulated "spokesmen" magnetic monopole 
introduced in the data (see sec. 2.1.1). 

For candidates with 4x 10~^ < (3 < 10^^ we compared 
the energy loss rate (computed with the assumption of 
a fixed 15 cm path length in the scintillator) with that 
expected for a MM (see Fig. ^). The /3 range was 

extended down to 4 x 10~^ to obtain an independent 
cross-check of the rejection based on the counter crossing 
time criterion for the candidates at the boundary between 
the two sub-intervals. In no case the energy deposition 
of a candidate was consistent with that expected for a 
monopole of the same apparent velocity. 

In this search magnetic monopole candidates were cat- 
egorized on the basis of the streamer tube and scintillator 
information. No event exhibits the pattern of a single par- 
ticle crossing the detector at the measured velocity; the 
velocity values were artifacts due to uncommon classes 
of events. A fraction of candidates 10 %) comes from 
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Ficenec et al. 




Fig. 4. Energy loss in the scintillator versus velocity for 
PHRASE events in the range 4 x 10"^ < P < 10"^ (for clarity 
only 40% of the candidates is shown). The calculated MM 
energy loss rates from [^,^ are also given, as a 95 % C.L. 
upper limit on the maximum energy loss of the candidates. 



occasional timing errors, which produced an apparent low 
value for the velocity of a cosmic ray muon crossing the ap- 
paratus. In a larger fraction of cases (~ 20%), the low ve- 
locity value comes from accidental coincidences, within the 
allowed time window, between radioactivity hits and/or a 
cosmic ray muon occurring in different scintillator layers. 
The remaining candidates are due to cosmic ray muons 
which crossed a scintillator layer and stopped just before 
reaching a second one. The decay electron produces a hit 
in this second layer, with a typical 2 ^s time delay (cor- 
responding to the muon life time) with respect to the first 
hit. 

This analysis was applied to the whole data set 
collected from October 1989 to June 1995 with the lower 
part of the detector, and with the full detector up to 
December 2000. The detector acceptance was computed 
by a Monte Carlo simulation, which takes into account 
the fraction of apparatus effectively in acquisition in each 
individual run. The total exposure was / Adt = 1.2 x 
10^^ cm^ssr. By taking into account the trigger and 
analysis efficiencies we obtain a 90 % C.L. flux upper limit 
<Z> < 2.2 X 10-16 cm-2 s-i sr"! for 1.2 x 10"^ < (3< lO-\ 
which is represented by curve "PHRASE" in Fig. K 



2.2 Searches with Streamer Tubes 

The streamer tube system was designed to be effective 
in the search for magnetic monopoles in a wide velocity 
range, 10^'' < (3 < 1. For this purpose we used a 
gas mixture containing helium and n-pentane. Helium is 
necessary to exploit the Drell effect the passage of a 
magnetic monopole with 1.1 x 10"'' < (3 < 10~^ through 
a helium atom leaves it in a excited metastable state. The 
n-pentane can then be ionized by collision with excited 
helium atoms (Penning effect). The high cross section for 



the processes ensures 100% efficiency in detecting bare 
MMs with this gas mixture. 

In the higher velocity region {(3 > 10^"^), where the 
assumptions used in the Drell-Penning effects do not 
apply, the standard ionization mechanism expected for a 
MM ensures an energy release several orders of magnitude 
higher than that due to minimum ionizing particles. The 
charge collected on a wire has a logarithmic dependence 
on the energy released inside the active cell 24 , so that 
a charge measurement allowed one to distinguish between 
MMs and muons. 

Horizontal streamer tubes were equipped with readout 
strips in order to provide a three dimensional event recon- 
struction. The memory depth of the readout electronics 
was large enough to store signals for monopoles with /3 
down to lO"'*. In order to retrieve spatial coordinates, sig- 
nals were shaped at ~ 550 /xs and sent to parallel- in/serial- 
out shift register chains for readout. Analog ORs of wire 
signals (covering 1 m for central planes and 0.5 m for lat- 
eral planes) were fed to the Charge and Time Processor 
(QTP) system for arrival time and charge measurement 
of the streamer signals . The trigger was based on the 
assumption that a heavy MM crossed the apparatus with- 
out any appreciable change in its direction and speed, due 
to its large kinetic energy. In designing the trigger logic, 
particular attention was given to avoiding the possibility 
that relativistic decay products from magnetic monopole 
induced nucleon decay, could cause dead-time in the ap- 
paratus before the MM generated the trigger ||25|| . 

Horizontal and vertical streamer tube planes were han- 
dled by two independent trigger systems. Monte Carlo 
simulations have shown that they were basically comple- 
mentary. Two analysis were performed using data col- 
lected with the two different triggers; all streamer tubes 
were used for the event reconstruction. Both analyses were 
based on the search for single tracks and on the measure- 
ment of the particle velocity by using the time informa- 
tion provided by the streamer tubes (maximum time jitter 
~ 600 ns) jlj]. Detailed investigations were performed in 
order to check that both trigger logics and analysis effi- 
ciencies were independent of the particle velocity. 



2.2.1 Search with horizontal streamer tubes 

This analysis technique is described in details in [Q. 
The assumpition made in the trigger design allows the 
trigger logic to be sensitive to any massive particle able to 
produce a signal in the streamer tubes. The streamer tubes 
time resolution allows identification of slow particles by 
measuring their time of flight across the apparatus ( "time 
track"). However, because of time jitter and afterpulsing 
in the streamer tubes, particles faster than 5 x lO^'^c could 
be confused with cosmic ray muons. 

The analysis was based on the search for single space 
tracks in both wire and strip views and on the measure- 
ment of the velocity of the candidates. The width of the 
temporal window is so large that on average three spuri- 
ous hits (due to the 40 Hz/m^ background on streamer 
tubes) are present in a 12 x 12 m^ plane per event. The 
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trigger logic selected events with an alignment in z (the 
vertical coordinate) versus time t of at least seven stre- 
amer tube planes in one or more of 320 time windows 
(called /3-slices). Each /3-slice covered a TOF window of 3 
yus and provided a rough estimate of the particle time of 
flight across the apparatus. 

The TOF value provided by the trigger allowed a 
discrimination of the streamer tube hits from a real track, 
from the spurious hits recorded in the ~ 500 /iS memory 
depth of the readout electronics. If a space track was 
found, a more refined alignment in the z — t temporal 
view was required ("time track"). This rejected effectively 
the background from radioactivity and the background 
from relativistic muons combined accidentally with some 
radioactivity hits. 

The muon provided the spatial track, while the ra- 
dioactivity hits may confuse the time tracking algorithm. 
Such events were rejected using a cut on the measured 
velocity. In fact, they were tagged with both the first (3- 
slice (that of relativistic muons) and with a higher /3-slice 
value (that of slow particles). For particles crossing the 
lower detector (path length ~ 5 m) the above condition 
corresponds to a cut for all particles with /3 > Pcut — 
(5 m/3 ns)/c = 5x10-3. 

The sensitivity of the trigger to relativistic muons 
provided a simple and efficient way to estimate the overall 
efficiency of this analysis. Muons crossing at least seven 
horizontal streamer tube planes survived all the analysis 
steps and were only rejected by the final /3-cut. The 
efficiency Sa of the analysis was then estimated by the ratio 
of the number of reconstructed muons to the expected 
number of single muons crossing at least seven planes. A 
further possible source of inefficiency were failures in the 
trigger circuits et ■ This effect was not accounted for by the 
above procedure, because muons were concentrated in the 
first /3-slice. The efficiencies £q and et, were completely 
described in [ p^ . 

This analysis used the data sample collected with the 
horizontal streamer tube trigger from January 1992 to 
September 2000, integrating a livetime of 8.1 years. The 
overall average efficiency was 74%. The detector accep- 
tance, computed by a Monte Carlo simulation which in- 
cluded geometrical and trigger requirements, was 4250 m^sr. 

No monopolc candidates were found. For 1.1 x 10^* < 
/3 < 5 X 10^3 the flux upper limit is 2.8 x 10" 
at 90% C.L. (Fig. |, curve "stream. H" 



2.2.2 Search with vertical streamer tubes 

Like in the search with horizontal streamer tubes, also in 
this case, the reconstruction procedure started from time 
and space information provided by the trigger. These were 
used to select hits compatible in both time and position 
with the /3-slice that fired. Using these time-position hits, 
we performed a complete space and time tracking of each 
events to reject accidental alignments; this leaves only 
the muon background. This was also used to evaluate 
the efficiencies. The spatial and temporal reconstruction 
algorithms are the same as those used in the horizontal 
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Fig. 5. The /? distribution for the candidate events recon- 
structed with the vertical streamer tube analysis. Events with 
P > Pcut = 3 x 10""^ are muons; all analyzed events are in- 
cluded. 



analysis. On the basis of the overall event reconstruction 
in the 3 spatial views (x — z, d—z, y—z) and in the {t — z) 
view, it was possible to achieve the complete geometrical 
and temporal reconstruction of the particle track and to 
compute its /3. The x — z spatial view corresponds to 
the horizontal streamer tube wire view; the d — z view 
corresponds to the horizontal streamer tube strips and 
the y — z view corresponds to that of the lateral (vertical) 
streamer tube wire view. Finally the t — z view is the 
temporal view along the vertical direction. As for the 
search with horizontal streamer tubes, muons are rejected 
if /3 > /3,"„t = 3 X 10-3 see Fig. | 

The analysis efficiency was estimated with the same 
approach used for the horizontal streamer monopole search. 

The search covered data from October 1994 to Septem- 
ber 2000 for a total of 4.4 liveyears. The overall average 
efficiency was 70%. The acceptance, estimated by Monte 
Carlo simulation, was 3018 m^sr. The measured (3 dis- 
tribution, shown in Fig. ^, is broader than the one ob- 
tained from the horizontal analysis (see ||l^ for compar- 
ison). This limited the sensitivity of this search to the 
velocity range 1.1 x 10^^ < /3 < 3 x lO^^. ]\[o magnetic 
monopole candidate was found. We thus establish an up- 
per limit to the monopole flux <P < 8x10^^^ cm^^ s^^ sr^-'^ 
for lO"-^ < /3 < 3 X IQ-^ (Fig. ||7curve "stream.V"). 



2.3 Searches with the Track-Etch detector 

The MACRO nuclear track subdetector covered 
a total area of 1263 m^ and was organized in stacks 
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Fig. 6. Measured reduced etch rate p vs REL for the CR39 
track etch detector exposed to relativistic and slow ion beams; 
the points are the experimental data, the solid line is the best 
fit to the data points {8 < Z < 68) obtained with relativistic 
ions. 



("wagons") of 24.5 x 24.5 cm^ consisting of three layers of 
CR39, three layers of Lexan and an aluminum absorber 
placed in an aluminized Mylar bag filled with dry air. 
The formation of an etchable track in a nuclear track 
detector is related to the Restricted Energy Loss (REL) 
Ip^ . There are two contributions to REL: the electronic 
energy loss (S'e), which represents the energy transferred 
to the electrons, and the nuclear energy loss {Sn), which 
represents the energy transferred to the nuclei in the 
material. In pTj it was shown that in our CR39 Sn is as 
effective as 5*6 in producing etchable tracks. The response 
of a nuclear track detector is measured by the reduced 
etch rate p = vt/vb, where vt and vb are the etching 
rates along the particle track and of the bulk material, 
respectively The reduced etch rate p vs REL for 

our CR39 is shown in Fig. it was measured using both 
relativistic and slow (down to /3 ~ 4 x 10"'^) ions [ p7| . 

The CR39 allowed a search for magnetic monopoles 
of different magnetic charges and /3. For a single Dirac 
charge, our CR39 is sensitive for MM in the ranges 
2 X 10-5 < /? < 2 X 10"'' and /3 > 1.2 x IQ-^ g. Lexan has 
a much higher threshold making it sensitive to relativistic 
MMs only (/3 > 10"^). 

Our CR39 was manufactured by the Intercast Europe 
Co. of Parma (Italy). A specific production line was set 
up in order to achieve a low detection threshold, high 
sensitivity in a large range of energy losses, high quality 
of the post-etched surface after prolonged etching and 
stability of the detector response over long periods of time 



We analyzed 845.5 ni^ of the MACRO CR39 track-etch 
sub-detector, with an average exposure time of 9.5 years. 
Since no candidates were found, the Lexan foils were not 
analyzed. The top CR39 foils were strongly etched in a 
8N NaOH water solution at 85° C till their final thickness 
reached 300-400 fim. The signal looked for was a hole or a 
biconical track with the two base cone areas equal within 
experimental uncertainties. 

After etching, the foils were scanned twice, using 
back lighting by different operators at low magnification 
looking for any possible optical inhomogeneity; the double 
scan guarantees an efficiency close to 100% for finding 
an etched hole. Detected inhomogeneities were further 
observed with a stereo microscope and were classified 
either as surface defects or as particle tracks. The latter 
were further observed under an optical microscope with 
high magnification. The axes of the base-cone ellipses in 
the top and bottom surfaces of the foils were measured and 
the corresponding p and incidence angle computed. A 
track was defined as a candidate if p and on the top 
and bottom sides were equal to within 15%. At a residual 
thickness of 300-400 /zm, double etch-pit tracks could be 
induced by proton recoils from neutron interactions or by 
low energy nuclei from muon interactions in the material 
surrounding the apparatus. 

For the few candidates remaining after the analysis 
of the first sheet (an average of 5/m^), we looked for 
a coincidence in position, angles and RELs in the third 
(bottom) CR39 layer. The second foil was etched in 6N 
NaOH water solution at 70 °C for 30 hours. An accurate 
scanning at high magnification was performed within an 
area of about 1 mm^ around the expected position of the 
candidate. 

We estimate a global efficiency of the procedure of 
99%. No two-fold coincidence was found, that is no 
magnetic monopole candidate was detected. The 90% 
C.L. fiux upper limits are at the level of 2.1 x 10""'^^ 
cm^^s^^sr^^ for g = go magnetic monopoles with /3 ^ 
10""* and 1.5x10"^^ cm^'^s^^sr^^ for magnetic monopoles 
with j3 ~ 1 (Fig. I curves "CR39"). 

2.4 Combined search for fast magnetic monopoles 

A search for fast magnetic monopoles with scintillators 
or streamer tubes alone was affected by the background 
due to high energy muons with large energy losses in the 
detector. This background could be reduced by a combi- 
nation of geometrical and energy cuts imposed on each of 
the sub-detectors. A combined analysis, based on all three 
subsystems, allowed the use of these cuts in a rather con- 
servative way. Moreover, any systematic error was greatly 
reduced by the combination of measurements from the 
three subsystems. The redundancy and complementary 
offered by the MACRO detector allowed a good rejection 
power against the background and a high reliability of 
possible candidates. The analysis procedure is fully de- 
scribed in jl^ . It used the data taken with the scintillator 
and the streamer tube sub-detectors to identify candidate 
events. This was done by reconstructing the energy release 
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Fig. 7. Combined searcli: Liglit yield for the seven MM 
candidates which survived the scintillator and streamer tube 
cuts. The candidates light yields are compared with the 
expected magnetic monopole signal as computed in (see 
Sect. 



(using the scintillators' Energy Reconstruction Processor 
-ERP- and the streamer tube Charge and Time Processor 
- QTP [HjH) and the particle's trajectory (using the strea- 
mer tubes' digital hit information). Any remaining events 
were then searched for in the track-etch layers as a final 
tool for their rejection or confirmation. The trigger selec- 
tion criteria required at least one fired scintillation counter 
and 7 hits in the horizontal streamer planes of the lower 
subdetector iQ. Once the event tracking was performed, 
the value of the energy lost in the scintillator intercepted 
by the track was reconstructed using the ERP system. The 
reconstructed energy in each selected scintillation counter 
was required to be AE > (Z\£')™„, = 150 MeV. The mi- 
nimum light yield by a MM in 10 cm pathlength in the 
scintillator is ~ 230 MeV. A further selection was applied 
on the streamer tube pulse charge by using the multiple 
measurements provided by QTPs along the particle tra- 
jectory. 

A cut |l4| was applied on the event average streamer 
charge, by exploiting the logarithmic dependence of the 
streamer charge on the primary ionization ||2^ , ^ . 

In the analyzed data set, seven events survived the 
cuts. For these events the track-etch wagons identified by 
the streamer tracking system were extracted and analyzed. 
No track compatible with the crossing of a MM was found. 
As a further check, for each candidate, the measured 
value of the energy lost in the scintillation counters 
was compared with the expected signal of a magnetic 
monopole of the same velocity [EsUlq]. This comparison 
is shown in Fig. 0. For three of the seven events, the 
TOF information was provided by the streamer system 
alone, since only one scintillation counter was present. 
In this case the error on the reconstructed velocity is 
large because of the limited time resolution (3 /is) of 
the streamer tubes. For all the events the measured 
energy losses were well below the expectations for a 



Fig. 8. Upper limits for an isotropic flux of magnetic mono- 
poles obtained by the different MACRO analyses as described 
in the text. The global limit is presented as curve "MACRO". 



MM. The analysis referred to about 4.8 live years; the 
maximum geometrical acceptance, computed by Monte 
Carlo methods, including all the analysis requirements, 
was 3565 m^ sr. The average global efficiency was 77%. 

Since no candidate survived, we set a 90% C.L. flux 
upper limit at 5.5 x 10~^^ cm~^ s~^ sr~^ (presented in 
Fig. ^ as curve "stream. -I- scint." ) for magnetic monopoles 
with 5 X 10-3 < P < 0.99. The analysis ends at /3 = 0.99 
since for Lorentz factor 7 10 magnetic monopoles might 
induce showers in the detector, reducing the analysis 
efficiency |14 . 



3 The global MACRO limit 

No magnetic monopole candidates were found in any 
of the above mentioned searches. Since any one of the 
sub-detectors may rule out, within its own acceptance 
and sensitivity, a potential candidate, a global limit was 
computed by combining together the limits obtained by 
the single analyses. The final MACRO limit and the 
limits obtained by the single subdetectors are given in 
Tab. |l| and shown in Fig. ^. Each search "j" produced 
a 90% C.L. flux limit given by <Pi = 2.3/Ai, where 
is the analysis time integrated acceptance. In order to 
obtain the global MACRO limit, we divided the full /3 
interval in a number of slices sufficient to characterize the 
changes in the individual acceptances. We required that 
the significance of the global limit (in terms of C.L.) is not 
altered. 

In order to illustrate the algorithm, suppose that in 
a specific (3 slice there are two analyses, "1" and "2", 
based on the use of two sub-detectors during the same 
period of time. We consider analysis "1" as "dominant" in 
the sense that it contributes with its full time integrated 
acceptance Ai to the global time integrated acceptance 
Aa- Analysis "2" then contributes to the global limit only 
with its independent parts relative to "1". We consider 
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Table 1. The 90% C.L. flux upper limits (in units of 10"'^'' cm"^ s"^ sr"^) as a function of /3 for an isotropic flux oi g = go 
magnetic monopoles with m > lO'^^ GeV/c^. The limits discussed in Section 2 are given in columns two to six; the global 
MACRO limit discussed in Section 3 is given in the last column. 



both the temporal independence as well as the spatial 
(geometric) independence versus the dominant analysis. 
The temporal independence is determined by comparing 
the "time efBciencies" e* of the analyses, defined as the 
ratios of each analyses live time to the covered solar 
time. If e\> the coefficient representing the temporal 
independence of "2" versus "1" is c| ^ = £3 — e*; otherwise, 
c\i = Q. In the case of the track-etch subdetector, there 
is no dead-time, so its temporal efficiency is equal to 1. 
The coefficient representing the geometric independence of 
analysis "2" versus "1", c| originates from the difference 
between the acceptances of the analyses. It is obtained by 
Monte Carlo simulations, assuming an incoming isotropic 
flux of magnetic monopoles with respect to subdetector 
"2": c^ i = {N2 - iVi,2)/A^2, where N2 and 7Vi,2 are 
the number of MMs detected by "2" and both analyses, 
respectively. 

The global time integrated acceptance is then: 
Ag = Ai+ 4 + (1 - 4,i)ci,i^2 

The global 90% C.L. limit for the flux of magnetic 
monopoles is = 2.3/Ag- 

The algorithm used to combine the actual MACRO 
limits is more complicated than the example above. 
For each analysis we took into consideration its actual 
history, eliminating the longer periods of time in which 
it was eventually missing, and the changes in the detector 
configuration (super- modules involved) . Those corrections 
were more critical in the case of earlier analyses, that 
were carried on during the construction of the MACRO 
detector and during initial tests; note that limits obtained 
by such older searches are not presented in Fig. ^ and in 
Tab.|l|, as they are considerably higher than the included 
ones, but they have still their imprint on the global limit. 

In Fig. I we present the global MACRO limit; for 
comparison, the flux limits from other experiments which 
searched for magnetic monopoles with similar properties, 
are also shown [ pl[ . In the figure the arrow indicates the 
Extended Parker Bound (EPB) at the level of 1.2 x IQ-^^ 
(m/10^^) cm~^s~^sr~"'^, which was obtained by consider- 



Parker Bound 



Ohya 




Fig. 9. The global MACRO limit for an isotropic flux of bare 
magnetic monopoles, with in > 10^'^ GeV/c^, g = go and 
o"cat < few mb. For comparison, we present also the flux limits 
from other experiments [Ell. 



ing the survival probability of a magnetic monopole of 
mass m in an early magnetic seed field p2| . 



4 Discussion 

Our analysis applies to an isotropic flux of bare MMs with 
charge g = gu=e/2a and nucleon decay catalysis cross 
sections smaller than 1 mb . 

The magnetic monopole flux at the detector site 
is isotropic if magnetic monopoles have enough kinetic 
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Table 2. Final results of magnetic monopole searches with 
the MACRO experiment. The limits depend on the magnetic 
monopole mass and on its velocity. 



5 Conclusions 

We present the final results of GUT magnetic monopole 
searches performed with the MACRO detector at Labo- 
ratori Nazionali del Gran Sasso (Italy). Different searches 
using the MACRO sub-detectors (i.e. scintillation coun- 
ters, limited streamer tubes and nuclear track detectors) 
both in stand alone and combined ways, were performed. 
Since no candidates were detected. Tab. 2 summarizes the 
90% C.L upper limit to an isotropic flux of bare MMs 
with charge g — g]j—e/2a and nucleon decay catalysis 
cross section smaller than 1 mb. 



energy (i.e. large mass and/or high /3) in order to cross 
the Earth. If this is not the case, only a fraction of the 
total detector acceptance is actually exploited and the 
upper limits given above must be corrected accordingly. 
In particular, if monopoles have sufficient energy to cross 
the overburden mountain (i.e. ~ 3700 hg/cm^) and reach 
the detector from above, but not enough to cross the 
Earth, the flux upper limit is about a factor of two 
weaker. Therefore the energy losses suffered by magnetic 
monopoles in traversing the Earth set the accessible values 
of mass and (3 for a given experiment. For MACRO, at 
least one half of the geometrical acceptance is ensured for 
relativistic MMs with m > 10^ GeV/c^ or, for /? ~ 10"" 
if m > 10^° GeV/c^; the full acceptance is reached if m > 
10^° GeV/c^ or m > 10"^^ GeV/c^, for relativistic and slow 
magnetic monopoles, respectively 34 (see Table 2). 



Sensitivity to fast and light MMs is also important, since 
they could be responsible for ultra high energy cosmic ray 
events above the Greisen-Zatsepin-Kuzmin cutoff 

If magnetic monopoles have g > go and/or an asso- 
ciated electric charge (dyons the detection efficiency 
might change. As far as our sub-detectors are concerned, 
the dyon detection efficiency would be greater than that 
of bare monopoles since the excitation/ ionization-based 
energy losses would be larger |^ . The only exception is the 
detection of very slow dyons with the streamer tubes, since 
it relies on the Drell effect on helium atoms |2^ , which for 
dyons might be prevented by coulombian repulsion. How- 
ever, as shown in ||^, the only effect would be the raising of 
the minimum velocity threshold at /? ~ (1.7 x 10~^)^/Z for 
dyons with electric charge Ze, the threshold for bare MMs 
being /3 = 1.1 x 10"'*. As suggested in [Q, GUT magnetic 
monopoles may catalyze nucleon decays along their path 
with a cross section tXc of the order of the hadronic cross 
sections [|l|. If CTc > 1 mb, the efficiencies of the aforemen- 
tioned searches might decrease due to the effects of the 
fast decay products. A deep study was performed on the 
effect of the nucleon catalysis cross section on the stre- 
amer tube analysis. The main result is that the present 
analyses are still efficient up to at least CTc — 100 mb. A 
dedicated search for magnetic monopoles accompanied by 
one or more nucleon decays along their path was also per- 
formed. The results are reported in a separate paper p7|. 
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